TOPOLOGICAL INSULATOR AND THE DIRAC EQUATION 
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We present a general description of topological insulators from the point of view of Dirac equations. 
The Z2 index for the Dirac equation is always zero, and thus the Dirac equation is topologically 
trivial. After the quadratic B term in momentum is introduced to correct the mass term m or the 
band gap of the Dirac equation, the Z2 index is modified as 1 for mB > and for mB < 0. For a 
fixed B there exists a topological quantum phase transition from a topologically trivial system to a 
non-trivial one system when the sign of mass m changes. A series of solutions near the boundary in 
the modified Dirac equation are obtained, which is characteristic of topological insulator. From the 
solutions of the bound states and the Z2 index we establish a relation between the Dirac equation 
and topological insulators. 

PACS numbers: 



INTRODUCTION 

Translational invariance in crystal lattices and the 
Bloch theorem for the wave function of electrons in solid 
make it possible for us to know the band structures of 
solid and why a solid is a metal, an insulator or semi- 
conductor. Recent years it is found that a new class 
of materials possess a feature that its bulk is insulating 
while its surface or edge is metallic. This metallic behav- 
ior is quite robust against impurities or interaction, and 
is protected by the intrinsic symmetry of the band struc- 
tures of electrons. The materials with this new feature is 
called topological insulator. [ll-Q 

In 1979 the one-parameter scaling theory predicted 
that all electrons in systems of two or lower dimension 
should be localized for a weak disorder. |4j This theory 
shaped the research of lower dimensional systems with 
disorders or interaction. Almost at the same time, von 
Klitzing et al discovered experimentally integer quantum 
Hall effect (IQHE) in two-dimensional (2D) electron gas 
in semiconductor hetero-j unction in a strong magnetic 
field, in which longitudinal conductance becomes zero 
while the quantum plateau of the Hall conductance ap- 
pears at ve 2 /h (y is an integer). [B[ Two years later Tsui 
et al observed the fractional quantum Hall effect (FQHE) 
in a sample with higher mobility. [f| In the theory of edge 
states for IQHE electrons form discrete Landau levels in 
a strong magnetic field. Electrons in the bulk has a van- 
ishing group velocity, and are easily localized by impu- 
rities or disorder while the electrons near the boundary 
are skipping along the boundary to form a conducting 
channel. [7J, |8| Thus in IQHE all bulk electrons are local- 
ized to be insulating while the edge electrons form several 
conducting channels according to the electron density 
which is robust against the impurities. This feature indi- 
cates that IQHE is a new state of quantum matter, i.e., 
one of topological insulators. In FQHE it is the electron- 
electron interaction that makes electrons incompressible 
and to form stable metallic edge states. 0, [l(| The quasi- 



particles in FQHE have fractionalized charges, and obeys 
new quantum statistics. In 1988 Haldane proposed that 
IQHE could be realized in a lattice system of spinless 
fermions in a period magnetic flux.[ll| Though the to- 
tal magnetic flux is zero, electrons are driven to form 
an conducting edge channel by the magnetic flux. Since 
there is no pure magnetic field the quantum Hall conduc- 
tance originates from the band structure of fermions in 
the lattice instead of the discrete Landau levels. 

In 2005 Kane and Mele generalized the Haldane's 
model to a lattice of spin 1/2 electrons. f]~2jj The strong 
spin-orbit coupling, an effect of relativistic quantum me- 
chanics for electrons in atoms, is introduced to replace 
the periodic magnetic flux in Haldane's model. This in- 
teraction looks like an spin-dependent magnetic field to 
employ on electron spins. Different electron spins ex- 
perience opposite spin-orbit force, i.e., spin transverse 
force. [l^| As a result, a bilayer Haldane model may be 
realized in a spin-1/2 electron system with spin-orbit cou- 
pling, which exhibits quantum spin Hall effect (QSHE). 
In the case there exist spin-dependent edge states around 
the boundary of the system: electrons with different spins 
move in opposite directions, and form a pair of helical 
edge states. In this system the time reversal symmetry 
is preserved, and the edge states are robust against im- 
purities or disorders because the electron backscattering 
in the two edge channel is prohibited due to the sym- 
metry. Bernevig, Hughes and Zhang predicted that the 
QSHE effect can be realized in the HgTe/CdTe sand- 
wiched quantum well. [3] HgTe is an inverted-band ma- 
terial, and CdTe is a normal band one. Tuning the 
thickness of HgTe layer may lead to the band inversion 
in the quantum well, which exhibits a topological phase 
transition. [15l 1 1 €31 ] This prediction was confirmed experi- 
mentally by Konig et al soon after the prediction. [TtJ The 
stability of the QSHE was studied by several groups. 0- 
EH ] Li et al discovered that the disorder may even gener- 
ate QSHE, and proposed a possible realization of topo- 
logical Anderson insulator, in which all bulk electrons are 
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localized by impurities meanwhile an conducting helical 
edge channels appear. 12211 T his phase was studied numer- 
ically and analytically. |23l l24j Strong Coulomb interac- 
tions may also generate QSHE in Mott insulators. 2|| 2(| 

The generalization of QSHE to three dimension is 
topologically non-trivial. 27-3(| Kane and Mele proposed 
a Z 2 index to classify the materials with time reversal in- 
variance into a strong and weak topological insulator. 31 
For a strong topological insulator, there exists an odd 
number of surface states crossing the Fermi surface of 
the system. The backscattering of electrons in the surface 
states are prohibited because of the symmetry. Bii_3;Sb x 
was predicted to be 3D topological insulator by Fu and 
Kane [32j and was verified experimentally. [33j| Zhang et al 
[3^ | and Xia et al[Ii| pointed out that Bi 2 Te3 and Bi 2 Se 3 
are topological insulator with a single Dirac cone of the 
surface states. ARPES data showed clearly the existence 
of single Dirac cone in Bi 2 Se 3 [35| and Bi 2 Te3 (36j. Elec- 
trons in the surface states possess a quantum spin texture 
structure, and electron momenta are coupled strongly 
with electron spins. These result in a lot of exotic mag- 
netoelectric properties. Qi et al. [37jproposed the uncon- 
ventional magneto-electric effect for the surface states, in 
which electric and magnetic fields are coupled together 
and are governed by so-call "axion equation" instead of 
Maxwell eqautions. It is regarded as one of the charac- 
teristic features of the topological insulators [H, [3^| . Fu 
and Kane proposed a possible realization of Majorana 
fcrmions as an proximity effect of s-wave superconduc- 
tor and surface states of topological insulator. [40} The 
Majorana fermions are topologically protected from lo- 
cal sources of decoherence, and will be of potential appli- 
cation in universal quantum computer. [4 ll |42| Thus the 
topological insulators open a new route to explore novel 
and exotic quantum particles in condensed matters. 

The Dirac equation is a relativistic quantum mechani- 
cal wave function for elementary spin 1/2 particle. 43ll44| 
It enters the field of topological insulator in two as- 
pects. First of all, topological insulators possess strong 
spin-orbit coupling, which is a consequence of the Dirac 
equation. [45| It makes the spin, momentum and the 
Coulomb interaction or external electric fields couple to- 
gether. As a result it is possible that the band structures 
in some materials becomes topologically non-trivial. An- 
other aspect is that the effective Hamiltonians to the 
QSHE and 3D topological insulators have the identical 
mathematical structure of the Dirac equation. In these 
effective models the equations are used to describe the 
coupling between electrons the conduction and valence 
bands, not the electron and positions in Dirac's the- 
ory. The positive and negative spectra are for the elec- 
trons and hole in semiconductors not in the high energy 
physics. The conventional Dirac equation is time-reversal 
invariant. For a system with time reversal symmetry, the 
effective Hamiltonian to describe the electrons near the 
Fermi level can be derived from the theory of invariants. 



As a result of the k ■ p expansion of the band structure, 
some effective continuous models have the same form of 
the Dirac equation. 

In this paper we start with the Dirac equation to pro- 
vide a simple but unified description for a large family 
of topological insulators. A scries of solvable differential 
equations are presented to demonstrate the existence of 
edge and surface states in topological insulators. 



DIRAC EQUATION AND SOLUTIONS OF THE 
BOUND STATES 



In 1928, Paul A. M. Dirac wrote down an equation for 
relativistic quantum mechanical wave functions, which 
describes elementary spin-1/2 particles, 43, 44| 
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cp ■ a 



mc j3 



(1) 



where m is the rest mass of particle and c is the speed of 
light, cti and (3 are the Dirac matrices satisfying 
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(2a) 
(2b) 
(2c) 



In 2D spatial space, the Dirac matrices have the same 
forms of the Pauli matrices cr.;, i.e., a x — a Xl a y = a y , 
and /3 = a z . In three dimensional spatial space, one 
representation of the Dirac matrices in terms of the Pauli 
matrices Oi (i = x, y, z) is 
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(3a) 
(3b) 



where ctq is a 2 x 2 identity matrix. From this equa- 
tion, the Einstein's relativistic energy-momentum rela- 
tion will be automatically the solution of the equation, 
E 2 = m 2 c A +p 2 c 2 . This equation demands the existence 
of antiparticle, i.e. particle with negative energy or mass, 
and predates the discovery of position, the antiparticle of 
the electron. It is one of the main achievements of mod- 
ern theoretical physics. 

Under the transformation of m — > —m it is found that 
the equation remains invariant if /3 — > —/3, which satisfies 
all mutual anticommutation relations for on and (3. This 
reflects the symmetry between the positive and negative 
energy particles. 

Possible relation between the Dirac equation and the 
topological insulator can be seen from a solution of the 
bound state at the interface between two regions of posi- 
tive and negative masses. For simplicity, we first consider 
a one-dimensional (ID) example 



h(x) — —ivhd x <j x + m(x)v 2 a z 



(4) 
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and 



m(x) — 



-mi if x < 
-m% otherwise 



(5) 



(and mi and 7772 > 0). Except for the extended solutions 
in the whole space, there exists a solution of the bound 
state with zero energy 



*(x) 



V 77117712 



\m(x)vx\/h 



(6) 



h mi + 7772 

The solution dominantly distributes near the point of 
x = and decays exponentially away from the point of 
x = 0. The solution of mi = 7772 was first obtained by 
Jackiw and Rebbi, and is a basis for the fractionalized 
charge in one-dimensional system.]^]. The solution ex- 
ists even when 7712 — > +00. In this case, Sfr(x) — > for 
x > 0. However, we have to point out that the wave func- 
tion does not vanish at the interface when 7712 — > +00. 
If we regard the vacuum as a system with an infinite 
positive mass, a system of a negative mass with an open 
boundary condition forms a bound state near the bound- 
ary. This is the source of some popular pictures for topo- 
logical insulator. 

In 2D, we consider a system with an interface parallel 
to the y-axis, with 777(2) = mi for x > 0, and —7772 for 
x < 0. k y is a good quantum number. We have two so- 
lutions which the wave functions dominantly distributes 
around the interface. One solution has the form 

a 



^(X,ky 



V 777 1 7772 
h mi + 7772 



-\m(x)vx\/h-\-ik y y 



(7) 



with the dispersion = vhk y 
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Another one has the form 



h mi 



777 2 



-\m{x)vx \ / h+ik y y 



(8) 
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with the dispersion = —vhk y . Both states carry a cur- 
rent along the interface, but electrons moving in opposite 
directions. The currents decays exponentially away from 
the interface. As the system does not break the time 
reversal symmetry, the two states are counterpart with 
time reversal symmetry with each other. This is a pair 
of helical edge (or bound) states at the interface. 

In 3D, we can also find a solution for the surface states, 
The dispersion relation for the surface states are e p = 
±vp. It has a rotational symmetry and forms a Dirac 
cone 

From these solutions we found that the edge states 
and surface states exist at the interface of systems with 
positive and negative masses. However, since there is a 
positive-negative mass symmetry in the Dirac equation, 
we cannot simply say which one is topologically trivial or 
non-trivial. Thus the Dirac equation alone is not enough 
to describe the topological insulators. 



MODIFIED DIRAC EQUATION AND Z2 
TOPOLOGICAL INVARIANT 

To explore the topological insulator, we start with a 
modified Dirac Hamiltonian by introducing a quadratic 
correction —Bp in momentum p to the band gap or 
rest-energy term, 



H 



vp ■ a - 



Bp 2 ) (3. 



(9) 



where mv 2 is the band gap of particle and m and v 
have dimensions of mass and speed, respectively, the 
quadratic term breaks the mass symmetry in the Dirac 
equation, and makes this equation topologically distinct 
from the original one. 

The general solutions of the wave functions can be ex- 
pressed as ^ u — u v (p)e l ( p ' r ~ Ep : l,t ' > / h . The dispersion rela- 
tions of four energy bands are E p i/ ( =12 ) = — ■Ep,i/(=3,4) — 
\Jv 2 p 2 + (mv 2 — Bp 2 ) 2 . The four-component spinors 
u v [p) can be expressed as u v [p) = Su v (p = 0) with 



2E, 



•p,i 



( 1 


p z v 

~ 

p+v 






p z v 


p—v 


£p 


Ep 


f 


P+v 


p z v 


e P 


e P 


P-V 


1 









p z v 





f 



(10) 
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where p± = p x ±ip y , t p = E, j 



p,i + (mi; 2 - Bp 2 ) , and u„(0) 
is one of the four eigen states of fi. 

The topological properties of the modified Dirac equa- 
tion can be gained from these solutions of a free particle. 
The Dirac equation is invariant under the time-reversal 
symmetry, and can be classified according to the Z2 topo- 
logical classification following Kane and Mele.[3lj]. I n 
the representation for the Dirac matrices in Eq. (??), 
the time- reversal operator here is defined as 47] & = 
—ia x a z K, where K. the complex conjugate operator that 
forms the complex conjugate of any coefficient that mul- 
tiplies a ket or wave function (and stands on the right 
of /C). Under the time reversal operation, the modified 
Dirac equation remains invariant, QH(p)Q~ 1 = H(—p) 
(p is a good quantum number of the momentum). Fur- 
thermore we have the relations that Qui(p) — —iu2(—p) 
and Qu2(p) — +iui(—p), which satisfy the relation of 
2 = — 1. Similarly, Qu^(p) — —iv,i(—p) and Qu^ip) = 
+iu^(—p). Thus the solutions of {ui(p), U2(—p)} and 
{u3(p), Ui(~ p)} are two degenerate Kramer pairs of pos- 
itive and negative energies, respectively. The matrix of 
overlap {(w M (p)| \u„(p))} has the form 
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■ mv 2 —Bp 2 




-Bp 2 





p z v 



P-V j p z v \ 
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• PzV 



; mv 2 — Bp 2 



E. 



pa 



&p a 
I p+ v 
Epa 

■ mv —Bp 2 
Bp a 
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which is antisymmetric, (u fi (p)\Q\u v (p)) = 
— {u v (p)\ For the two negative energy 

bands u^(p) and u^(p), the submatrix of overlap can be 
expressed in terms of a single number as ea V P(p), 



P(p) = i 



Bp 2 



\J (mv 2 — Bp 2 ) 2 + v 2 p 2 



(12) 



which is the Pfaffian for the 2x2 matrix. According to 
Kane and Mele,[3l| the even or odd number of the ze- 
ros in -P(p) defines the Z2 topological invariant. Here we 
want to emphasize that the sign of a dimensionless pa- 
rameter mB will determine the Z2 invariant of the mod- 
ified Dirac equation. Since -P(p) is always non-zero for 
mB < and there exists no zero in the Pfaffian, we con- 
clude immediately that the modified Dirac Hamiltonian 
for mB < including the conventional Dirac Hamilto- 
nian (B = 0) is topologically trivial. 

For mB > the case is different. In this continuous 
model, the Brillouin zone becomes infinite. At p — and 
p = +00, P(0) = isgn(m) and P(+oo) = — zsgn(P). In 
this case P(p) = at p 2 = mv 2 /B. p = is always 
one of the time reversal invariant momenta (TRIM). As 
a result of an isotropic model in the momentum space, 
we can think all points of p — +00 shrink into one point 
if we regard the continuous model as a limit of the lattice 
model by taking the lattice space a — > and the recip- 
rocal lattice vector G = 2n/a —> +00. In this sense 
as a limit of a square lattice other three TRIM have 
P(0,G/2) = P(G/2,0) = P(G/2,G/2) = P(+oo) which 
has an opposite sign of P(0) if mB > 0. Similarly for 
a cubic lattice P(p) of other seven TRIM have opposite 
sign of P(0). Following Fu, Kane and Mele[27], we 
conclude that the modified Dirac Hamiltonian is topolog- 
ically non-trivial only if mB > 0. 

In two dimension Z2 index can be determined by eval- 
uating the winding number of the phase of P(p) around 
a loop of enclosing the half the Brilouin zone in the com- 
plex plane of p = p x + ip y , 



i = ^-.i d P -Vpiog[P(p) + iJ]. 

2m J c 



(13) 



Because the model is isotropic, the integral then reduces 
to only the path along p x -axis while the part of the half- 
circle integral vanishes for 5 > and |p| — > +00. Along 
the p x axis one one of a pair of zeros in the ring is enclosed 
in the contour C when mB > 0, which give a Z2 index 
1=1. This defines the non-trivial QSH phase. 



TOPOLOGICAL INVARIANTS AND QUANTUM 
PHASE TRANSITION 



An alternative approach to explore the topological 
property of the Dirac model is the Green function 
method, (iil Volovik [ijj proposed that the Green func- 
tion rather than the Hamiltonian is more applicable to 



classify the topological insulator. From the Dirac equa- 
tion, the Green function has the form 

G{iuJ n ,p) 



iui n — H 
vp ■ a + (mv 2 — Bp 2 )(3 — ius n 

where h 2 (k) = H 2 = v 2 p 2 + (mv 2 — Bp 2 ) 2 . There is the 
following topological invariant 

N = -^ei jk Tr[Ki Un = dpGd Pi G- 1 Gd Pj G- 1 Gd Pk G- 1 } 

where K = a y <g> a$ is the symmetry-related operator. 
After tedious algebra, it is found that 

N = sgn(m) + sgn(P). 

When mB > 0, N — ±2, which define the phase topo- 
logically non-trivial. If B is fixed to be positive, there 
exist a quantum phase transition from topologically triv- 
ial phase of m < to a topologically non-trivial phase. 
This is in a good agreement with the result of Z2 index 
in the preceding section. 

Except for the phases of N = ±2, it is found that 
there exist a marginal topological phases of JV = ±1. 
For free Dirac fermions of B = 0, the topological invari- 
ant N = sgn(m). It is +1 for a positive mass and -1 
for a negative mass. Their difference AN = 2 which is 
the origin of the existence of the bound states at the in- 
terface of two systems with positive and negative mass 
as we discussed in Section II. There exists an interme- 
diate gapless phases of m = between two topological 
nontrivial (N = ±2) and trivial (N = 0) phases. At the 
critical point of topological quantum phase transition, all 
intermediate states are gapless. Its topological invariant 
is also N = +1 or — 1 just like as the free Dirac fermions. 



THE TOPOLOGICALLY PROTECTED 
BOUNDARY STATES SOLUTIONS 

ID: the bound state of zero energy 

Let us start with the ID case. In this case, the equation 
in Eq. ([9]) can be decoupled into two sets of independent 
equations in the form 



h(x) = vp x a x + (mv 2 - Bp 2 ) a z 



(14) 



For a semi-infinite chain, we consider an open boundary 
condition at x = 0. We may have a series of extended 
solutions which spread in the whole space. In this section 
we focus on the solutions of bound states near the end 
of the chain. We require that the wave function vanishes 
at x = +00. In the condition of mB > 0, there exists a 
solution of the bound state with zero energy 



^\ = C_( sgn(P) 
V2 V 1 



) (15) 
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id 



2\B\h 



(i±VT 



4mi? 



(16) 



where C is the normalization constant. The main feature 
of this solution is that the wave function dominantly dis- 
tributes near the boundary. The two parameters £_ > £ + 
and decides the spatial distribution of the wave function. 
This is a very important length scale, which character- 
izes the bound state. When B —> 0, £ + — > \B\ H/v and 
£_ = h/mv i.e., £ + approaches to zero, and £_ becomes 
a finite constant. If we relax the constraint of the van- 
ishing wave function at the boundary, the solution exists 
even if B = 0. In this way, we go back the conventional 
Dirac equation. In this sense, the two equations reach at 
the same conclusion. 

In the four-component form to Eq.©, two degenerate 
solutions have the form, 



C_ 
3 



/ sgn(S) \ 



) (17a) 



*2 = 



c_ 

V2 



( \ 

sgn(S) 

i 

V o J 



( e -*/«+_ e -*/«-) (17b) 



2D: the helical edge states 

In two dimension, the equation is decoupled into two 
independent equations 



vp x a x ± vpydy + (mv 2 



Bp 2 ) a z 



(18) 



These two subsets of equations breaks the " time" reversal 
symmetry under the transformation of oi — > —<Ji and 
Pi -> ~Pi- 

We consider a semi-infinite plane with the boundary 
at x = 0. p y is a good quantum number. At p y = 0, 
the 2D equation has the same form as the ID equation. 
The x-dependent part of the solutions of bound states 
has the identical form as in ID. Thus we use the two ID 
solutions {^ r i,4 r 2} a s the basis. The y-dependent part 
Ai?2D = vpyOty — Bpyfl is regarded as the perturbation to 
the ID Hamiltonian. In this way, we have a ID effective 
model for the helical edge states 



H eff = ((* 1 \ > (* 2 \)AH 



l*i> 



vp y sgii(B)a z (19) 



The sign dependence of B in the effective model also 
reflects the fact that the helical edge states disappear if 
B = 0. The dispersion relations for the bound states at 
the boundary are 



Electrons will have positive (+v) and negative velocity 
(— v) in two different states, respectively, and form form 
a pair of helical edge states. Thus the 2D equation can 
describe a quantum spin Hall system. 

The exact solutions of the edge states to this 2D equa- 
tion have the similar form of ID 1 50] 



C_ 

V2 



*2 = —!= 



c_ 
7i 



( sgn(B) \ 



i 



sgn(B) 



-*/«, 



with the dispersion relation E Py = ±vp y sgn(B)a z . The 
penetration depth becomes p y dependent, 



€±1 = 2\Bjh i 1 ± y/l-*™B + 4B*p*/v< 



(22) 



In two-dimension, the Chern number or Thouless- 
Kohmoto-Nightingale-Nijs integer can be used to charac- 
terize whether the system is topologically trivial or non- 
trivial. [H| Write the Hamiltonian in Eq. (|T8"1) in the form 
H = d(p) ■ a The Chern number is expressed as 



, F .. u( ]. diLdfc 
dp tjjfc^i g Px q Pv 



d 3 



where d 2 = J2 a =x y z ^q-0> [52|| The integral runs over 
the first Brillouin zone for a lattice system. The num- 
ber is always an integer for an finite first Brillouin zone, 
but can be fractional for an infinite zone. For these two 



equations the Chern number has the form [53|, [54 1 



n± = ±(sgn(m) + sgn(£?))/2. 



(23) 



which gives the Hall conductance a± — n±e 2 /h. When m 
and B have the same sign, n± becomes ±1, and the sys- 
tems are topologically non-trivial. But if m and B have 
different signs, n± = 0. The topologically non-trivial 
condition is in agreement with the existence condition of 
edge state solution. This reflects the bulk-edge relation 
of integer quantum Hall effect. (HBI 



3D: the surface states 

In 3D, we consider an y-z plane at x = 0. We can 
derive an effective model for the surface states by means 
of the ID solutions of the bound states. Consider p y - and 
p z -dependent part as a perturbation to ID Hid(x), 



±VPy 



(20) 



AH; 



3D 



vpydy + vp z a z - B(p 2 + p 2 z )j3. (24) 
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The solutions of 3D Dirac equation at p y — p z = are 
identical to the two ID solutions. A straightforward cal- 
culation as in the 2D case gives 

H eff = , (*,\)AH 3D ( ) = v S gn(B)(p x a) x . 

(25) 



Under a unitary transformation 
1 



<I>, - -^(l*!)-^)) 



(26a) 



$ 2 = -= (|*i) +i|* 2 » (26b) 
V2 



APPLICATION TO REAL SYSTEMS 

Now we address the relevance of the modified Dirac 
model to real materials. Of course we cannot simply ap- 
ply the Dirac equation to semiconductors explicitly. Usu- 
ally the band structures of most semiconductors or others 
have no particle-hole symmetry. Thus a quadratic term 
should be introduced into the modified Dirac model. On 
the other hand the band structure may not be isotropic 
and the effective velocities along different axes are differ- 
ent. A more general model has the form, 

H = efa) + ViPiVi + ( mv2 ~ B ^)P- ( 3 °) 



we can have a gapless Dirac equation for the surface 
states 



H, 



eff = i((1>l|,($ 2 |)A^(j^ 

= vsgn(B) (p y a y +p z a z ). 



(27) 



The dispersion relations become E p — ±vp. In this way 
we have an effective model for a single Dirac cone of the 
surface states. 

The exact solutions of the surface states to this 3D 
equation with the boundary are 

*± - CV° ± (e- x /t+ - e-*lt- ) exp[+i ( Py y + p z z) /ft] 

(28a) 

where 



To have a solution for topological insulator, the addi- 
tional terms must keep the band gap open. Otherwise it 
cannot describe an insulator. 

The equation in solids can be derived from the theory 
of invariant or the k-p theory as an expansion of the mo- 
mentum p near the T point. Since under the time rever- 
sal, P — ¥ j3 and a — > —a, if we expand an time reversal 
invariant Hamiltonian near the T point, the zero-order 
term should be constant, e(0) and mv 2 f3. The first order 
term in the momentum must ViPiOti since pi — > — pi 
under time reversal. The second order term is Bip^/3 

2 

and w— in e(pi). The third order term is the cubic term 
in a. The summation up to the second order terms give 
the modified Dirac equation. 



Complex p-wave spinless superconductor 



/ cos |sgn(_B) \ 
—i sin |sgn(B) 



2 k 

sin ■ 
i cos ■ 



\ i cos | / 

( sinfsgn(B) ^ 
icos|sgn(B) 
— cos J 
\ i sin 5 ,/ 



(28b) 



(28c) 



with the dispersion relation E± = ±upsgn(_B) and p = 
p 2 + p 2 z . The penetration depth becomes p dependent, 



2\B\h 



1 ± y/l-4mB + 4B 2 p 2 /fr 



(29) 



A complex p-wave spinless superconductor has two 
topologically distinct phases, one is the str ong pa iring 
phase and another is the weak pairing phase. 48|,|56| The 



weak pairing phase is identical to the Moore-Read quan- 
tum Hall state. [56( The system can be described by the 
modified Dirac model. In the BCS mean field theory, the 
effective Hamiltonian for quasiparticles in this system has 
the form 



K, 



eff 



E 

k 



^4 c fc + \{K c -kCk + A fe 4d fc ) 



The normalized ground state has the form 
\n)=H(u k+ v k clcl k )\0). 



(31) 



(32) 



Generalization to higher dimensional topological 
insulators 



The solution can be generalized to higher dimensional 
system. We conclude that there always exists a d- 
dimensional surface state in the modified Dirac equation. 



where |0) is the vacuum state. The Bogoliubov-de 
Genncs equation for u k and v k becomes 



ih 



d 



Uk 



dt V v k 



-Afe 



(33) 



For complex p-wave pairing, we take to be an eigen 
function of rotations in k of eigenvalue of two-dimensional 
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angular momentum! I = — 1, and thus at small k it gener- 
ically takes the form 



A k = A(k w -ik y );Z k = — -11 (34) 



2m 



In this way the Bogoliubov-de Gennes equation has the 
exact form of 2D modified Dirac equation 



H e ff = -A {k x a x + kytjy) + ( — - fi ) a z (35) 



The Chern number of the effective Hamiltonian becomes 



n = [sgn(M) + sgn(l/m)] /2 



(36) 



Since we assume the mass of the spinless particles m pos- 
itive, we conclude that for a positive /i(> 0) the Chern 
number is +1 and for a negative \i the Chern number is 
0. For fx = 0, the Chern number is equal to one half, 
which is similar to the case of m —> +oo and a finite \i. 
If the quadratic term in £ fc is neglected, we see that the 
topological property will change completely. 

Usually for a positive //, the system is in a weak pair- 
ing phase, for a negative fi the strong coupling phase. 
Including the quadratic term in £ fc we conclude that the 
weak pairing phase for positive /i is a typical topologi- 
cal insulator, which possesses a chiral edge state if the 
system has a boundary. The exact solution of this equa- 
tion can be found in the paper by Zhou et al.[5(| Read 
and Green argued that a bound state solution ex- 
ists at a straight domain wall parallel to the y-axis, with 
fx(r) = fx(x) small and positive for a; > 0, and negative 
for x < 0. There is only one solution for each k v and 
so we have a chiral Majorana fermions on the domain 
wall. From the 2D solution, the system in a weak pairing 
phase should have a topologically protected and chiral 
edge state of Majorana fermion. Recently Fu and Kane 
proposed that as a superconducting proximity effect the 
interface of the surface state of three-dimensional topo- 
logical insulator and an s-wave superconductor resembles 
a spinless p x + ip y superconductor, but does not break 
time reversal symmetry, (ioj The state support Majorana 
bound states at vortices. 



Quantum Spin Hall Effect: HgTe/CdTe quantum 
well and thin film of topological insulator 

In 1988 Haldane proposed a spinless fermion model 
for IQHE without Landau levels, in which two indepen- 
dent effective Hamiltonian with the same form of 2D the 



Dirac equation were obtained. [Ill] The Haldane's model 
was generalized to the graphene lattice model of spin 1[2 



electrons, which exhibits quantum spin Hall effect. [12 1 
Bernevig, Hughes and Zhang predicted that QSHE can 
be realized in HgTe/CdTe quantum well and proposed 



an effective model, 14 1 
Hbhz 



_ ( h{k) 
~ I h*(-k) 



(37) 



where h(k) = e(k) + A(k x a x + k y a y ) + (M - Bk 2 )a z . 
The model is actually equivalent to the 2D Dirac model 
as shown in Eq. (fT8)) in addition of the kinetic term e(fc), 

h(k) = e(k) + h+; h*{-k) = e(k) + Uh-U' 1 , (38) 

where the unitary transformation matrix U = a z . 

If the inclusion of e{k) does not close the energy gap 
caused by M for a non-zero B, there exists a topologi- 
cal phase transition from a positive M to a negative M. 
However, the sign of M alone cannot determine whether 
the system is topologically trivial or non-trivial. From 
the formula in Eq. (|23[) . we know that the system is in 
a quantum spin Hall phase only for MB > and there 
exists a pair of helical edge states around the boundary 
of system. A general discussion can be found in the pa- 
per by Zhou et al.jHoj Finally we want to comment on 
one popular opinion that the band inversion induces the 
topological quantum phase transition. If B = 0, the sys- 
tem is always topologically trivial for either positive or 
negative M, though there exists a bound state at the 
interface of two systems with positive and negative M, 
respectively. 

The surface states of a thin film of topological insu- 
lator such as Bi 2 Te3 and Bi 2 Se 3 can be also described 
by a two-dimensional Dirac model. 5^, 54 1 The mass or 
the band gap of the Dirac particles originates from the 
overlapping of the wave functions of the top and bottom 
surface states. The gap opening of the two surface states 
were observed in Bi2Se3 thin films experimentallv|57l |58| . 
and was also confirmed numerically by DFT[59i]. Re- 
cently Luo and Zunger [6(| reported a DFT calculation 
for HgTe/CdTe quantum well and presented a different 
picture that the topological quantum phase transition 
occurs at the crossing point of two " interface- localized" 
states. This is in a good agreement of the theory for 3D 
topological insulator thin film. 61 1 



Three-Dimensional Topological Insulators 

The 3D Dirac equation can be applied to describe a 
large family of three-dimensional topological insulators. 
Bi2Te3 and Bi2Se3 and Sb2Te3 have been confirmed to 
be topological insulator with a single Dirac come of sur- 
face states. For example, in Bi2Te3, the electrons near 
the Fermi surfaces mainly come from the p-orbitals of Bi 
and Te atoms. According to the point group symmetry 
of the crystal lattice, p z orbital splits from p XjV orbital. 
Near the Fermi surface the energy levels turn out to be 
the p z orbital. The four orbitals are used to construct 
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the eigenstates of parity and the base for the effective 
Hamiltonian, 1 34] which has the exact form as 



H = e(k) + \ ViPiOti + (mv 2 



i=x 1 y,z 



E ( 39 ) 



i—x,y,z 



with v x 



v n and v z 



v± and B T 



By = _B|, 



and B z = B±. e(k) = C - D u {p 2 x + p 2 ) - D±p 2 . In this 
way the effective Hamiltonian in the x-y plane has the 
form 53] 



Heff = y/1 - Dl/Blv„(p x a) 2 
or under a unitary transformation 



(40) 



H eff = yjl - D 2 jB\v H {p x a x + Py a y ). (41) 

We note that the inclusion of e(k) will revise the effective 
velocity of the surface states, which is different from the 
result in Ref . [34 1 . 



FROM THE CONTINUOUS MODEL TO THE 
LATTICE MODEL 

In practice, the continuous model is sometimes mapped 
into a lattice model in the tight binding approximation. 
In a d-dimensional hyper-cubic lattice, one replaces [6^, 
H 



hi 



sin kj a 



a 
2 



(1 — COS fcjflj) 



(42) 
(43) 



CONCLUSION 

To summarize, we found that the Z2 index for the Dirac 
equation is always zero, and thus the Dirac equation is 
topologically trivial. After the quadratic B term is intro- 
duced to correct the mass m of the Dirac equation, the Z2 
index is modified as 1 for mB > and for mB < 0. For 
a fixed B there exists a topological quantum phase tran- 
sition from a topologically trivial system to a non-trivial 
one system when the sign of mass m changes. 

From the solutions of the modified Dirac equation, we 
found that under the condition of mB > 0, 

• in ID, there exists the bound state of zero energy 
near the boundary; 

• in 2D, there exists the solution of helical edge states 
near the boundary; 

• in 3D, there exists the solution of the surface states 
near the surface; 

• in higher dimension, there always exists the solu- 
tion of higher dimension surface. 

From the solutions of the bound states near the 
boundary, and the calculation of Z2 index we conclude 
that the modified Dirac equation can provide a descrip- 
tion of a large families of topological insulators from one 
to higher dimension. 
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which are equal to each other in a long wave limit. Usu- 
ally there exits the fermion doubling problem in the lat- 
tice model for massless Dirac particles. The replacement 
of ki — > sin ki a/ a will cause an additional zero point at 
fcja = 7r besides kia = 0. Thus there exist two Dirac 
cones in a square lattice at k = (0,0) and (Tr/a,w/a) 
for a gapless Dirac equation. A large B term removes 
the problem as 2B(1 — cos kia)/a 2 — > B/a 2 in the lattice 
model. Thus the lattice model is equivalent to the con- 
tinuous model only in the condition of a large B. The 
zero point of (1 — cos^a) 2 is at fcj<z = ir/2 not or ir 
in sinfcja. Thus for a finite B, the band gap may not 
open at the T point in the lattice model because of the 
competition between the linear term and the quadratic 
term of ki . This fact may lead to a topological transition 
from a large B to a small B. Imura et al [gij analyzed 
the 2D case in details and found that there exists a topo- 
logical transition at a finite value of B in two-dimension. 
A similar transition will also exist in higher dimension. 
It should be careful when we study the continuous model 
in a tight binding approximation. 
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